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Abstract 

The linear and non-linear evolution of a relativistic current sheet of pair (e^) plasmas is inves- 
tigated by three-dimensional particle-in-cell simulations. In a Harris configuration, it is obtained 
that the magnetic energy is fast dissipated by the relativistic drift kink instability (RDKI). How- 
ever, when a current-aligned magnetic field (the so-called "guide field") is introduced, the RDKI 
is stabilized by the magnetic tension force and it separates into two obliquely-propagating modes, 
which we call the relativistic drift-kink-tearing instability (RDKTI). These two waves deform the 
current sheet so that they trigger relativistic magnetic reconnection at a crossover thinning point. 
Since relativistic reconnection produces a lot of non-thermal particles, the guide field is of critical 
importance to study the energetics of a relativistic current sheet. 

PACS numbers: 52.27.Ep, 52.27.Ny, 52.35.Vd, 95.30.Qd 



It is widely believed that plasma heating and particle acceleration occur in a wide variety 
of plasma regions that contain magnetic fields. A current sheet structure where the magnetic 
field polarity changes its direction is one of the most fundamental structures among them. 
Importantly, when two groups of magnetic field lines with opposite polarities meet each 
other around the current sheet, magnetic reconnection takes place and then it explosively 
dissipates the magnetic energy into the kinetic energy of plasmas. In fact, reconnection is 
accepted as a main player in stellar and solar flares [1] and storms in the Earth's magne- 
tosphere [2J- The theory of reconnection has often been studied in a current sheet with 
anti-parallel field lines but it is possible that field lines are somehow "twisted" and so the 
system with a finite amplitude of a current-aligned magnetic field (the so-called "guide 
field" ) has recently been investigated by three-dimensional (3D) simulations [31 HI El El [7] . 
Magnetic reconnection processes are also important in high-energy astrophysical contexts; 
the jets from active galactic nuclei [5J [91 [TP] , pulsar wind [HI [12] and probably gamma-ray 
bursts [131 E] . Particularly in the Crab Nebula, it has been a long-standing problem (the 
so-called "a-problem" ) how originally Poynting-dominated plasmas [15J are converted into 
kinetic-dominated plasmas [16] in the relativistic outflow from the neutron star. The most 
promising solution is relativistic magnetic reconnection [17] in the striped current sheets 

However, reconnection processes in such a relativistic hot e ± plasmas (the plasma tem- 
perature T is larger than the rest mass energy: T 3> mc 2 ) has been poorly understood. 
Recently, 2D particle simulations of relativistic magnetic reconnection in such a relativistic 
current sheet (hereafter RCS) were carried out [TBI HE]- O ne important point is that a lot of 
magnetic energy can be converted into non-thermal energy of plasmas due to the enhanced 
dc acceleration around the X-type region. By the way, a RCS is extremely unstable to the 
relativistic drift kink instability (RDKI) [20] , a relativistic extension of the drift kink insta- 
bility (DKI) [2U 1221 123] , whose wave vector is normal to 2D reconnection plane. Contrary 
to the relativistic reconnection, the RDKI mainly dissipates magnetic energy into plasma 
heat rather than non-thermal energy. In addition, the effect of the guide field in a RCS has 
never been investigated in 3D/2D. 

We have performed 3D particle-in-cell simulations to explore linear /nonlinear evolution 
of a RCS, firstly taking into account the guide field effect. The system consists of 256 3 
grids. Periodic boundaries in the x, y and z directions and with double current sheets 
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are considered in the z direction. The half width of the current sheet A is set to 10 grids 
so that the boundaries are located at x — ±12.8, y = ±12.8 and z = ±6.4 in unit of 
A. Also, time is normalized by the light transit time r c = A/c. We take a relativistic 
Harris model as the initial current sheet configuration [12] . The magnetic field and their 
distribution functions are described by B x = £> tanh(;z), B y = aB Q , f± oc n(z) exp[—T ^{e — 
(3±mcu y }/T] + ribg exp(— e/O.lT) and n(z) = n cosh~ 2 (z) , where B is the magnitude of 
anti-parallel component of the magnetic field in the lobe (the background region), a is the 
relative amplitude of the guide field to the reversed field, n is the number density of plasmas 
in the current sheet, n^g is the number density of plasmas in the lobe, f3± = v±/c are the drift 
velocities for each species; (3+ = +(3 for positrons, (3- = —(3 for electrons, Tp = [1 — /3 2 ] -1 / 2 , 
e is the particle energy and u is the four velocity of u = [1 — (v / c) 2 ] -1 / 2 ■ v. We investigate 
two cases of a = (Run A) and a = —0.5 (Run B). We set T = mc 2 and f3 = 0.3. The total 
number of super particles is 5 x 10 s and we set the plasma density to Uq ~ 80 (pairs) and 
n bg ~ 4 (pairs) per one grid so that nb g /n = 0.05. The total energy is conserved within an 
error of 0.5% throughout the simulation runs. 

First, we present a result of Run A without guide field. A snapshot at t = 80 is presented 
in panel (a) in Fig. [l] The system evolution is similar to that of the 2D RDKI [20] , including 
profiles of perturbed physical properties. The wave-number of the most dominant mode is 
k y ~ 0.74 (mode 3) and its growth rare is u>i/Q c = 0.02, where Q c is the gyro-frequency; 
Q c = uj c / r y = (eB ) j ('-/mc) , 7 is the Lorenz factor for typical particle energy. In the 2D 
work, the growth rate has its peak around k y ~ 0.7 — 0.8 at Ui/fl c = 0.03 — 0.04. After 
this stage, the RCS becomes folded, turns into the non-linear stage and then mixed into 
the turbulent state. Because of the limited simulation height in z, some plasmas are also 
mixed across the z boundaries. Finally, more than 83% of the magnetic energy is dissipated 
into the particle energy. Fig. [2] shows energy spectra of 3D runs. For Run A, three stages 
(t = 0,80, 140) are presented. In later time of t = 140, one can recognize a high-energy tail 
due to the particle acceleration inside folded current sheet. More importantly, plasmas are 
heated throughout the late-time mixing stage of the RDKI. The transportation across the 
z boundaries provides no significant change between 2D and 3D in the spectra. We think 
that the observed growth rate becomes slower than that of 2D due to the limited system 
size in z. Note that it is still faster than that of 2D reconnection (oji/Vt c = 0.011). Since the 
RDKI grows faster when T > mc 2 [2U], we conclude that the RDKI is a dominant process 
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in a RCS. Again, magnetic dissipation and plasma heating by the RDKI would be the main 
signature of a RCS. 

Next, we show results of Run B with a finite amplitude of the guide field B y = —0.5B . 
Snapshots at two stages (t = 120 and 170) are presented in panels (b), (c) and (d) in Fig. [TJ 
Instead of the RDKI, a flute-like mode whose wavevector is k\ = (0.25,0.25), is observed 
on the upper side of the RCS at t = 120 in Fig. [I] (b). Similar flute mode is observed on 
the lower side, too, but its wavevector is in different direction: k 2 = (0.25, —0.25). In order 
to take a clear look at this lower-side mode, we present cross-sections at z = — 1 in panels 
(b) and (c). Snapshots at t — 170 are presented in panels (c) and (d) in Fig. [T] Due to 
the two flute modes, the RCS is so modulated that a plasma density hole appears at the 
center. In panel (d), cross-sections of the RCS at x,y,z = and typical magnetic field 
lines are presented. Magnetic reconnection occurs around the center of the simulation box. 
The speed of reconnection jets is up to 0.6c and the outflow plasmas are transported into 
the O-point(s) around (x,y) ~ (±12.8,0). Then, since magnetic field mainly consists of the 
guide field component B y inside the RCS, they are dissipated along the O-line region around 
x = ±12.8. The dense region around (x,y) ~ (0, ±12.8) is the remnant of the thick point 
of the RCS, where two obliquely-propagating modes are linked. Energy spectra in Run B 
are also presented in Fig. [2} The oblique mode itself has few effect on the spectra, however, 
after relativistic reconnection starts, strong enhancement of the non-thermal tail is observed 
due to the particle acceleration around the X-type region [TBI US]- 

Along with 3D work, we have studied an effect of the guide field using several sets of 
2D particle simulations (|a| = 0,0.25,0.5,1.0). It is obtained that reconnection in the xz 
plane grows slower when larger |a| is set, while we were unable to observe the the RDKI in 
the yz plane when |a| exceeds some critical value (|ot| > 0.5). We found that this is due 
to the magnetic tension effect. When a = 0, magnetic field lines are always parallel to the 
wavefronts of the instability and so they cannot be bent (no tension force) but when the 
guide field is introduced, the RDKI will be opposed by the tension force of current-aligned 
field lines. 

In order to study the obliquely-propagating modes in Run B, we investigate eigen func- 
tions and their linear growth rate of instabilities in the RCS. The relativistic two-fluid 
equations are linearized, assuming that perturbations are given by 5f oc f(z) exp(ik x x + 
ik y y — iut). The obtained growth rate (ui/Q, c ) are presented in contour maps in Fig. [3] for 
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two cases of a = 0, —0.5. In these maps, reconnection (tearing mode) is plotted along the k x 
axis and the RDKI is plotted along the k y axis. We mention that all of the obtained modes 
are purely growing and that growth rate for (k x ,k y ) are equal to the rate for (\k x \, \k y \) 
because of the symmetry. In Fig. [3] (a), the RDKI or its neighbor of the relativistic drift 
sausage instability (RDSI) has the maximum growth rate in k x = while we obtain the 
dominant mode at (k x , k y ) ~ (0, 0.74) in the 3D simulation. We note that these instabilities 
are dumped by the kinetic effect in shorter wavelength of \k\ > 1, where the fluid theory is 
not valid. Because of the similarity of the eigen profiles and growth rates, it is natural to 
say that the oblique modes are intermediate mode between the RDKI/RDSI and the tear- 
ing instability. Hereafter we shall call them the "relativistic drift-kink-tearing instability" 
(RDKTI) [21], the oblique extension of the RDKI/RDSI driven by fc-aligned component 
of the current. Since the RDKTI mode is weakly stabilized by the field lines and since 
its driving force is weaker, the RDKTI grows slower than relevant RDKI/RDSI in a RCS 
with exact anti-parallel fields. When we set larger |a|, the RDKI/RDSI along the k y axis 
become slower due to the magnetic tension of stronger guide field. Instead, two branches 
of the RDKTI for k = (k x ,±k y ) become dominant, or, in other word, the RDKI/RDSI 
separates into two branches of oblique RDKTI waves. Fig. [3] (b) shows growth rate in case 
of a = —0.5. As obtained by supplemental 2D runs, the RDKI along the k y axis is stabilized 
by the magnetic tension. The right panel in Fig. [4] shows eigen profiles of density perturba- 
tion for A; = (0.25, ±0.25). They are highly asymmetric in z; a mode for ki = (0.25,0.25) 
has its peak in the upper side of the RCS while the other one for fc 2 = (0.25, —0.25) has 
its peak in the lower side. Such asymmetry is explained by "twisted" effect of background 
magnetic fields. In case of the mode for k\, the magnetic field is quasi-perpendicular to fci 
in the upper side of the RCS so that its tension can be neglected. On the contrary, the field 
is quasi-parallel to k x in the lower side, so that the mode is opposed by strong magnetic 
tension. The RDKTI for k 2 is an upside-down mirror of the RDKTI for k\ because of the 
physical symmetry. As a result, we observe two dominant RDKTI modes: one for k\ in the 
upper side and the other for k 2 in the lower side of the RCS. As larger \a\ is set, an angle 
between two oblique wave- vectors (such as fci and k 2 ) becomes wider in accordance with 
the lobe magnetic field lines. In such cases, the RDKTI grows slower because of stronger 
magnetic tension. 

The secondary relativistic reconnection is triggered by the RDKTI in its non-linear stage. 
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Since two major RDKTI modes are in different direction, the deformed RCS become very 
thin at the cross-over point where both the upper and the lower RDKTI wave press the 
RCS. This mechanism is illustrated in Fig. |4j In addition to this enhanced thinning effect, 
accompanying electric fields E y and the heated plasma due to the RDKTI leads to the trig- 
gering of magnetic reconnection. Once reconnection broke up, it continue to grow and then 
particle acceleration takes place near the reconnecting region. It produces more nonthermal 
particles due to the direct particle acceleration in the reconnection region [T8J EES]. Actu- 
ally there are two crossover thinning points - one is around (x,y) ~ (0,0) and the other is 
around (x, y) ~ (±12.8, ±12.8). However, reconnection jets from the former point are finally 
transported into the latter point so that the reconnection does not dominate at the latter 
point. The speed of the reconnection jets is finally up to 0.78c at later stage. 

The electromagnetic energy (Eem) is dissipated by the rate of ^Eem ~ — ^VtLxLyB^/ Sir, 
where L x>y are the system size and v t is a typical dissipation speed [121 03] • We obtain 
v t ~ 0.25c in Run A and this value is consistent with the z-displacement speed of the RDKI: 
v z = 0.2 — 0.3c [20] • In Run B, v t ~ 0.05c or larger. The ratio of the total field energy 
{Eem) to the total particle energy (Ep = ^jmc 2 ) is originally 1.0 but it falls below 0.1 at 
t = 160 in Run A, while the ratio evolves from 1.25 to 0.7 at t = 220 in Run B. Although 
this ratio extremely depends on the system-size, it seems that the RDKI is favorable to 
produce kinetic-dominated plasmas. 

Finally, we discuss our system size limitations. The most dominant RDKTI has mode 
(1, ±1) in Run B. Since smaller-wavelength modes will be suppressed by the kinetic effect and 
since longer- wavelength RDKTI modes has smaller growth rate in Fig. [3] (b), the obtained 
modes would be roughly correct and as long as a pair of the RDKTI modes dominate, 
the trigger mechanism of reconnection should be the same. The global evolution of the 
RCS still remains unclear. One possibility is that multiple reconnection regions evolve 
into one big reconnection region, and then higher energetic particles will be produced in 
the wider acceleration region. Another possibility is that many reconnection regions grow 
and that the system becomes turbulent, which may contribute to statistical acceleration or 
thermalization. It should be further investigated by larger particle simulations or two-fluid 
relativistic MHD simulations. 

Let us summarize this paper. Considering the effect of the guide field, we present the 
following two scenarios in a RCS. One is plasma heating by the RDKI in an anti-parallel case 
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without the guide field. The other is non-thermal acceleration by secondary reconnection in 
a twisted case with the guide field. In this case, magnetic reconnection is triggered by the 
coupling of two RDKTI waves, which is separated from the RDKI/RDSI. It should be noted 
that slight difference in magnetic field topology have great influence on the destination of 
the magnetic energy: plasma thermal energy or non-thermal energy. 
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FIG. 1: (color) (a) Snapshot of the current sheet in Run A (a = 0) at t = 80. Two gray surfaces 
show n = 2/2>riQ. The plasma density at the neutral plane (z = 0) is projected into the bottom roof, 
with color from black (empty) to red (dense; n ~ 1.2no). The white arrow shows plasma flow. The 
light speed (v = c) is projected to the length of 4. (b) Snapshot of the current sheet in Run B 
with a guide field configuration (a = —0.5) at t = 120. Two gray surfaces show n = 2/3no and 
the plasma density under the neutral plane {z = — 1) is presented in color in the bottom roof, (c) 
Snapshot at t = 170 in Run B. (d) Another snapshot at t = 170 in Run B, with typical magnetic 
field lines. Three projections show cross-sections at x,y, z = 0. 
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FIG. 2: Energy spectra observed in Runs A and B. The initial state (the same in both two runs) 
and two typical stages for each runs are selected. Time is normalized by r c and particle energy is 
normalized by mc 2 . 
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(a) a =0.0 



(b) a =-0.5 




FIG. 3: (color online) Growth rate (u;j/Q c ) of the instabilities in wave-vector spaces of (0 < k x < 
1,0 < ky < 1). (a) a = and (b) a = -0.5. 
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Triggering of magnetic reconnection 




FIG. 4: Left: A schematic illustration of obliquely-propagating modes and the triggering mech- 
anism of magnetic reconnection. Right: eigen profiles of the relativistic drift-kink-tearing modes. 
A gray line stands for k\ = (0.25,0.25) and a light gray line shows its partner: k<i = (0.25, —0.25). 
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